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The study of electron-positron collisions at center-of-mass energies corre- 
sponding to charmonium production has reached new levels of sensitivity 
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Q_i! thanks to experiments by the BES and CLEO Collaborations. Final states 

\ jP, where P is a pseudoscalar meson such as it , 77, and 77' can arise ei- 

I ther from charmonium decays or in the continuum through a virtual photon: 

e + e~ — > 7* — > 7P. Estimates of this latter process are given at center-of- 
mass (cm.) energies corresponding to the J/i/j(1S), if)(2S), and t/>(3770) 
resonances. 
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gn| PACS numbers: 13.40. Gp, 13.66.Bc, 14.40.Aq, 14.40.Cs 



I. INTRODUCTION 



The decays of charmonium states to a photon and a neutral pseudoscalar meson 
P = (7T°, 77,77') can shed light on a number of different mechanisms, including two-gluon 



couplings to qq states, vector meson dominance, mixing of heavy quarkonium with bound 
states of light quarks and antiquarks, and final-state radiation by light quarks [Ij. In 
crj ■ these final states there can be direct contributions from the continuum process 

O 

g: e + e" -»• 7* -> 7P , (1) 

where P is a neutral pseudoscalar meson. The 7*-7-P vertex was shown [2] to be 
characterized by a form factor F(Q 2 ), where Q 2 = —q 2 and q is the four-momentum 
of the virtual photon 7*, behaving as F(Q 2 ) — > 2f 7T /Q 2 . Here f n = 93 MeV is the 
neutral-pion decay constant. In the present article we evaluate the cross sections for 
P = 7T°, 77, and 7/ at y/s = 3097, 3696, and 3773 MeV. These results are relevant to 
possible continuum backgrounds in searches for decays of J/ip, ^(3686), and ?/>(3770) to 
7P, currently undertaken by the CLEO Collaboration pp. 

In Section II we calculate the cross section for the process 01) , given the form factor 
estimate in Ref. [2]. We then estimate bounds on the actual cross section in Section III 
by comparing the expected result for e + e~ — > J/ip —>■ 777° using vector dominance with 
the observed (slightly larger) cross section. We summarize in Section IV. 

II. CONTINUUM 7P PRODUCTION 

The 7*-7-P vertex was discussed in Ref. [2J, where a proposal was made to test 
it in photon-photon collisions with one photon highly off-shell. In that case, the off- 
shell photon will have a spacelike Q 2 > 0. Experimental tests in this regime were indeed 
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performed in Ref . [3] . However, the reaction (CEJ) is equally useful in testing that vertex for 
timelike Q 2 = —q 2 < 0. The BaBar Collaboration [4] has measured the cross sections for 
e + e~ — > 7//') at q 2 = 112, GeV 2 providing valuable information on the j*^y-r]^ vertex 
in the asymptotic regime. We shall compare our predictions with their results presently 
The vertex of interest has the form 

= -te 2 F 7T o(Q 2 )^u P ,p u e p q C7 . (2) 

Here the form factor F n o{Q 2 ) is expected to behave for large Q 2 as [2] 

F n o(Q 2 )^2f n /Q 2 . (3) 

The four-momentum of the pseudoscalar meson (here, 7r°) is denoted as p, while e p is 
the polarization vector of the outgoing on-shell photon. 

The differential cross section with respect to cos#, where 9 is the angle the outgoing 
photon makes with the beam axis in the e + e~ center-of-mass system, is 

da{e+e - ,. ^ ^J^V Kp = 1 _MlL i (4) 

d(cosO) \ s J s 

where s = q 2 is the square of the center-of-mass energy and Mp is the mass of the 
pseudoscalar meson (here, 7r°). Integrating with respect to cos 8, we find 



a(e + e -> 7 -> 77r u ) = — ) K* . (5) 



This may be compared with the cross section for muon pair production (neglecting 



( + - * + -\ 47rc * 2 86 - 8 nb fa\ 

° {e e ^ 7 ^ " ) = ^ = JjGeV 2 -y (6) 
which is 9.05 nb at yfs = 3.097 GeV. Thus 

a(e+e — > 7* — ► ) s 

which is 4.11 x 10~ 5 at y/s = 3.097 GeV. At this energy we thus have 

a(e + e~ -> 7* -> 7tt°) = 372 fb . (8) 

The ratio (J7J) can also be used to predict the virtual-photon contribution to the 
branching fraction B(J/ip —>■ 77?°) in terms of B(J/ip — ► fi + fi + ): 

With B(J/if> — > n + assumed to be equal to B(J/ip — > e + e~) = 5.94% [5], this implies 
B{ J/ip -> 7* -> 7vr°) = 2.44 x 10" 6 , far below the observed value [5] of (3.3lg;|) x 10" 5 . 
We shall see in the next section that the vector-dominance process J /if} — > p°*7r° — > 77r° 
accounts for most if not all of the observed branching fraction. 



2 



Whereas the asymptotic form factor F 7T o(Q 2 ) = 2f n /Q 2 of Ref. [2] implies R^j/^ = 
4.11 x 10~ 5 , a slightly larger value of R&rj/M = 10~ 4 is implied by the calculation of Ref. 
[B|. This would imply 

a(e + e~ -> 7* -> 7tt°) = 905 fb . (10) 

for the continuum contribution at ^/s = 3.097 GeV. 

We shall now scale the cross section estimates (jHJ) and (TIP]) to other energies and 
particle types. Let the quark model wave function of a neutral particle P be characterized 
by the sum of pairs q^ with coefficients cf : 

l^> = E^I««>- (ii) 

i 

Then the ratio of the 7*-7-P form factor to that for ir° is just 

f p (q 2 )/faq 2 ) = Y,c(QVY,cfQi ■ (i 2 ) 

i i 

The 77 and 7/ may be represented as octet-singlet mixtures, 

77 = cos ^r/g + sin Qr\\ , 7/ = — sin #t) 8 + cos 0r/i , (13) 

r/ 8 = (uu + dJ— 2ss)/v^6 , rji = (uu + dd + ss) / V3 . (14) 

An approximate form which fits the data well [7], and which we shall take in what 
follows, is 

77 = (uu + dd — ss)/V3 , 7/' = (uu + dd+ 2ss)/VQ , (15) 

corresponding to a mixing angle 8 = — arcsin(l/3) = —19.5°. This implies 
\F v (Q 2 )/F n o(Q 2 )\ 2 = 32/27, |F v (Q 2 )/F 7r o(Q 2 )| 2 = 49/32. (For comparison, the unmixed 
octet and singlet states give \F m (Q 2 ) / ' F^(Q 2 )\ 2 = 1/3, |^ 1 (g 2 )/F 7r o(g 2 )| 2 = 8/3.) 
The above mixing ansatz implies 

linv^lg 5 '-F V (Q 2 )\ = 2fJ - = 202 MeV , (16) 

Hm g2 ^ 00 |g 2 F ?7 ,(g 2 )| = 2/^- = 251 MeV, (17) 

in satisfactory agreement with the results from BaBar [4 J at q 2 = 112 GeV 2 : 

| g 2 F r ,(Q 2 ) I = 229 ± 30 ± 8 MeV , (18) 
|g 2 F^(Q 2 )| = 251 ±19 ±8 MeV. (19) 

Other ansdtze for 77-77' mixing have been explored, for example, in Ref. [8j. The two-gluon 
components of 77,77' are not considered here, but are treated in Ref. [9j. 

The kinematic factors K 3 are summarized for P = 7r°, 77,77' and cm. energies corre- 
sponding to J/if), ip(2S), and 7/^(3770) in Table fl] Also shown are the continuum cross 
sections at these cm. energies. At 3.77 GeV, the predicted continuum cross sections for 
e + e - _^ 77? (0 are (0.19,0.25) pb, consistent with the values of (0.17l^±0.03, 0.21±g;g|d= 
0.03) pb observed by CLEO [TJ. Thus, it is consistent to assume that the 777 and 777' 
signals at 3.77 GeV come entirely from continuum. Upper bounds on £>[t/>(3770) — > 777W] 
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Table I: Kinematic suppression factors K 3 = [1 — (Mj>/s)] 3 and continuum cross sections 
cr(e + e~ — > 7* — > 7P) for neutral pseudoscalar mesons P. 







K 3 






(fb) 




p 




if)(2S) 


^(3770) 


J/if; 


if)(2S) 


^(3770) 


7T° 


0.994 


0.996 


0.996 


372 


186 


169 




0.909 


0.935 


0.938 


403 


207 


189 




0.740 


0.811 


0.819 


502 


274 


252 



Table II: Predicted ratios Rg and 7P branching ratios for quarkonium states decaying 
to 7P via a virtual photon 7*. 



p 


j/lf) 


K (10- 

iP(2S) 


- 5 ) 

^(3770) 


J/if) 


B(jP) 
iP(2S) 


^(3770) 




4.11 


2.91 


2.78 


2.44 x 10~ 6 


2.19 x 10~ 7 


2.7 x 10~ 10 


V 


4.45 


3.24 


3.10 


2.64 x 10" 6 


2.43 x 10~ 7 


3.0 x 10~ 10 


T}' 


4.68 


3.62 


3.49 


2.78 x 10" 6 


2.73 x 10" 7 


3.4 x 10~ 10 



under various scenarios of interference between direct decay and interference are quoted 
in Ref. p. 

The ratios Rg characterize the branching ratios of charmonium states to 7P via 
the virtual photon 7*. We summarize these ratios and the corresponding predicted 
contributions to quarkonium branching ratios in Table [Til Here we have taken [5] 
B(J/ip -> n+fi-) = 5.94 x 10" 2 , B(ifj(2S) -»■ = 7.52 x 10~ 3 , and B(ip(377Q) -»■ 

fi + H~ = 9.71 x 10~ 7 . Except for the case of J/if) — > 77?°, which we shall discuss further 
in the next section, these contributions are negligible. 

All the estimates presented above were for the form factor F n o(Q 2 ) behaving as 
2f n /Q 2 [2]. The calculation of Ref. [6] would give Pg = 10 -4 at the J /if), and hence all 
values approximately 2.43 times as large. 

III. COMPARISON WITH DATA AND VECTOR DOMINANCE 

The observed branching ratios [5] for J/ if) and if)(2S) decaying to 7P are summarized 
in Table IIHI The Particle Data Group averages [5] include determinations by the BES 
Collaboration at the J /if) [TU] but not the recent CLEO results [1J. As mentioned 
earlier, the single- virtual-photon process J /if) — > 7* —>■ 77?° is insufficient to account for 
the observed branching ratio. 

We now show that most of not all of the observed B(J/if) — > 771°) can be accounted 
for by the vector dominance model (VDM). We follow the calculation of Ref. [6], in 
which 

T(J/lf) -» 77T )y D M = / P^o\ ( efffi 

T(J/iP^pH°) - [p^oj \m p 
Here the cm. 3-momenta in J/if) —>■ 77r° and J /if) — > p°7r° are p* = 1546 MeV/c 2 and 
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Table III: Observed branching ratios for [J/ip, ip(2S)] — > 7P [U [3] , in units of 10 4 . 



p 






^(25) 






Ref. [TJ 


Ref. 


Ref. [TJ 


Ref. [5J 


7T° 


0.363 ±0.036 ±0.013 


n qq+U.Uti 
u -°°-0.04 


< 0.07 


< 54 


V 


11.01 ±0.29 ±0.22 


9.8 ± 1.0 


< 0.02 


< 0.9 


rf 


52.4 ±1.2 ±1.1 


47.1 ±2.7 


1.19 ±0.08 ±0.03 


1.36 ±0.24 



1448 MeV/c 2 , respectively. The neutral p meson decay constant is obtained from that 
of the charged p meson using r decays. Specifically (see., e.g., Ref. [TT] ) 



fp± = V^/tt 



B(t — > v T p 



ml — m 2 m T 



m 2 _ m 2 



P \l m r ± 2 ™p 



(21) 



jB(t — > z/ r Ti 

= (209 ± 1.6) MeV (22) 
where particle masses and branching fractions are taken from [5]. Then 

f p0 = f p ±/V2 = 147.8 ± 1.1 MeV . (23) 
A very similar result is obtained from the expression 

T(P° - e+e-) = ^ /p 2 o . (24) 

based on the matrix element (0| J^ m |p°(g, e)) = e M m p / p o. With £>(p° — > e + e~) = (4.71 ± 
0.5) x 10~ 5 , m p = (775.49 ± 0.34) MeV/c 2 , and T p = (149.4 ± 1.0) MeV, one finds 
f p0 = (156 ± 8) MeV. 

We choose not to use the Particle Data Group average for B(J/ip — > pir), but to 
average the three most recent determinations of this branching fraction, summarized in 
Table HVl Taking 1/3 of this average for the p°7r° final state, we have B{J/ip — > p°7r°) = 
(7.10 ± 0.03) x 10~ 3 . Eq. (TJOJ then implies 

B(J/ip -> 7^°)vdm = (28.8 ± 1.3) x 10" 6 . (25) 

This value is consistent with the observed branching fraction quoted in Table 1111} but 
also allows for an additional contribution. 

In Ref. [B], the relative signs of the VDM and virtual-photon (7*) contributions are 
predicted to be positive, so that we may place rather restricted upper bounds on the 
latter. We take a 90% confidence level (c.l.) bound of B(J/vp — > 771- ) < 41 x 10~ 6 . 
The prediction of Ref. [6] was that S 7 *(J/-0 — ► r y7i°)/B(J/tp — > e + e~) = 10~ 4 , or (using 
B{J/ip -> e+e") = 5.94%), 

B r (J/ij) -> 7vr ) cz = 5.94 x 10- 6 (26) 

Let A be the maximum allowed fraction of the CZ branching ratio. Then the CZ assump- 
tion of constructive interference between the VDM and virtual-photon contributions 
implies 

(V28.8± 1.3 + V5.94 A) 2 < 41 , (27) 
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Table IV: Values of B(J/ip — ► pn) used in computing average, in percent. 



Source 


Reference 


Value 


BES direct J/ip 


m\ 


2.184 ±0.005 ±0.201 


BES J/ip from ip(2S) 


[12j 


2.091 ±0.021 ±0.116 


BaBar radiative return 


[13] 


2.18 ±0.19 


Average 




2.13 ±0.09 



or, taking the lower limit of the VDM theoretical error, A < 0.23. The CZ continuum 
cross section for e + e~ —>■'-/*—>■ 771- at the J/ip cm. energy, corresponding to 10~ 4 of 
that for muon pair production, is 905 fb, so the assumption of constructive interference 
between VDM and virtual photon contributions in J/ip — > 771" implies 

a(e + e~ -> 7* -> 7vr°) < 210 fb , (28) 

a bit more than half the value predicted in Table [J from the form factor in Ref . [2] . The 
agreement between prediction and data for 777W production at 3.77 GeV suggests that 
not all of the predictions of Table I are subject to the same suppression. 

It is quite possible that as \q 2 \ — > 00, the asymptotic values of the relevant form 
factors are approached from below. This is indeed what is found for spacelike q 2 < in 
singly-tagged photon-photon collisions where one photon is highly virtual |3j . Moreover, 
the corrections to the form factors in perturbative QCD are of the form 1 — [5a s (q 2 )}/ (3ir) 
[T4] , easily entailing a suppression of the cross section by 25% or more. 

IV. CONCLUSIONS 

We have presented continuum cross sections for e + e~ — > 7* — > 7P, where P is a 
neutral pseudoscalar meson 7r°, rj, or r]'. Calculations are presented at cm. energies 
corresponding to the masses of J/ip, ip(2S), and ^(3770) for the 7*-7-7r° form factor 
advocated by Brodsky and Lepage |2], and rescaled to the 77 and rf using their anticipated 
quark content and known kinematic factors. 

It is shown that the contributions of the 7*-7~P vertex to the branching fractions 
for the above-mentioned charmonium states to 7P are negligible except in the case of 
J/ip — ► 77T°, where the assumption of constructive interference with the VDM contribu- 
tion permits rather stringent bounds to be placed, equivalent to continuum cross sections 
for e + e~ — > 7* — ► 7P a bit more than half those based on the asymptotic form factor 
behavior predicted in Ref. [2]. 

The cleanliness of the CLEO-c detector environment has permitted new studies of 
7P final states at the J/ip, ip(2S), and ^(3770), with 777W signals at 3.77 GeV consistent 
with continuum production. 
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